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Exchange interaction within a quantum dot strongly affects the transport through it in the Kondo 
regime. In a striking difference with the results of the conventional model, where this interaction 
is neglected, here the temperature and magnetic field dependence of the conductance may become 
non-monotonic: its initial increase follows by a drop when temperature and magnetic field are 
lowered. 
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The conventional theory of the Kondo effect in tunnel- 
ing |l| is based on the Anderson impurity model Q . The 
use of this model for describing the electron transport 
through a conducting grain was first suggested in Q , and 
later successfully applied in the context of quantum dots 
in the Coulomb blockade regime (see, e.g., 0). In these 
applications, the quantum dot is modeled by one singly 
occupied energy level. The model predicts monotonic 
increase of the tunneling conductance with lowering the 
temperature. By now, a strong experimental evidence 
exists for such behavior 

To measure transport across the dot, at least two leads 
should be attached to it. The problem of transmission 
through a dot is similar to that of transition between 
channels in a multichannel scattering problem. If the 
system is described by the Anderson model, then the cor- 
responding transition rate turns out to be proportional 
to the total scattering cross-section. Therefore the two 
problems, calculation of the tunneling conductance G, 
and evaluation of the magnetic impurity contribution to 
the resistivity of a bulk metal, are equivalent ||. Well be- 
low the Kondo temperature, the scattering cross-section 
reaches the unitary limit; accordingly, G saturates at the 
level G ~ e 2 /nh. 

In a real quantum dot, the truncation of the spectrum 
to a single level is not always possible. For example, dots 
with ground state spin S > 1/2 are not described by the 
Anderson impurity model [Q. In this paper we demon- 
strate that in general the problem of transport through 
a dot does not map onto the problem of resistivity in a 
bulk metal. In the most interesting case, the tunneling 
conductance first rises and than drops when temperature 
is lowered. The conductance dependence on the magnetic 
field is also non-monotonic. 

The confining potential forming a lateral quantum dot 
is smooth. The dot-lead junction is essentially an elec- 
tron waveguide [gj. Making it narrower, one pinches the 
propagating modes off. Coulomb blockade develops when 
the last propagating mode is near its pinch-off. There- 
fore, for lateral dots in the Coulomb blockade regime the 
number of channels per junction is one. 

We start with a discussion of the dependence of the 
zero-temperature conductance on the magnetic field ap- 
plied in the plane of the dot. The in-plane field results 
in the Zeeman splitting B of the spin states of the dot, 



but barely affects the orbital degrees of freedom. At a 
finite B, the ground state of the system is not degen- 
erate; therefore, at T = an electron experiences only 
potential scattering. The amplitudes of scattering S s - aa > 
of electrons with spin s from lead a to lead a' form the 
scattering matrix (here a, a' = R or L for the right or 
left lead, respectively). It can be diagonalized by rota- 
tion in the R — L space to the new basis of channels a 
and b, 

US S U^ = diag {e 2i5lB } , j = a,b, s = ±1. (1) 

Here U = exp(idT y ) exp(i(pT z ), and t 1 are the Pauli ma- 
trices acting in the R— L space (t + = t x +ir y transforms 
L — > R). The parameters $ and (p depend on the micro- 
scopic properties of the dot-lead junctions. The T = 
conductance is given by the Landauer formula, 
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and can be expressed now in terms of the scattering phase 



shifts (5 7S at the Fermi level: 
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(note that the conductance is independent of <p) . 

In order to calculate the phase shifts, one needs to 
know the effective Hamiltonian Tt of the system. Clearly, 
the term in it representing the interaction of itinerant 
electrons with the spin of the dot should obey the SU(2) 
symmetry, 
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where s 7 = J2kk> ss > i J \ks( eT sa> /^)ipjk>s' are local spin den- 
sities of itinerant electrons. The last term in TL describes 
the effect of the magnetic field. The Hamiltonian (^) 
does not include the potential scattering term; we defer 
the discussion of its role till later in the paper. The de- 
scription by Eq. (^) is valid at energies below the energy 
gap A for spin excitations in the dot, which is of the or- 
der of the single-particle level spacing in it. Note that the 
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symmetry allows for any signs of the exchange constants 
J 7 . We will present a microscopic derivation of Eq. (^) 
towards the end of the paper. 

Following P], we remove the ambiguity in the defini- 
tion of the phase shifts by fixing <L, S = at J 7 = 0; 
then \5~f S \ < n/2 at any J 7 . Furthermore, the invari- 
ance of the Hamiltonian (||) with respect to the particle- 
hole transformation, ip.yks —> sipl / _ k _ s , yields the relation 
iS s iS_ s = 1 for the scattering matrix, which allows us to 
represent the phase shifts as 



Sj S = — s<5 7 , < (5 7 < 7r/2. 
Substitution into Eq. (g) then gives 
G/G Q = sin 2 (5 a ~S b ). 
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In order to find 8y at B = 0, we need to invoke the 
properties of the ground state of the Kondo model. We 
start with S > 1/2 on the dot. In the case of an anti- 
ferromagnetic exchange in a channel, J 7 > 0, the itiner- 
ant electrons participate in the screening of the localized 
spin, thus reducing its value by 1/2. By the Friedel sum 
rule, the corresponding phase shift (5 7 = 7r/2. The chan- 
nels with the ferromagnetic coupling to the localized spin, 
J 7 < 0, decouple at low energies, and <5 7 = 0. It is then 
clear from Eq. (||) that for S > 1/2 the conductance at 
T — and B = is enhanced due to the Kondo effect 
only if the two exchange constants have opposite signs. 

At S = 1/2 the enhancement of the conductance oc- 
curs also if both exchange constants are positive. Indeed, 
the local spin is screened by the channel, say a, with the 
larger exchange constant, so that <5 a = tt/2. The other 
channel decouples |l(| at low energies, and Sb = 0. Note 
that in the special case J a — <h the conductance across 
the dot is zero, which is obvious from the rotational sym- 
metry in the R — L space, existing in this case. 

There is no mapping between the Kondo problems for 
tunneling and for scattering in the bulk, except if one of 
the exchange constants in Eq. (|]) is zero. However, if the 
constants are finite and have opposite signs (or both are 
positive at S = 1/2), the conductance behavior is qual- 
itatively similar to the conventional || one: G{T) and 
G(B) increase monotonously when T and B are lowered, 
eventually saturating at the value G ~ e 2 /irfr. If both 
coupling constants are negative (fully ferromagnetic ex- 
change interaction) there is no Kondo effect: G decreases 
when T and B are lowered. 

We concentrate now on the most interesting case of 
the antiferromagnetic coupling, J a > Jb > 0, and spin 
S > 1. In this case, the Kondo effect manifests itself in 
a peculiar way: the field and temperature dependence of 
G is not monotonic, which differs qualitatively form the 
conventional picture. Here we present a detailed analysis 
only of the dependence G(B) at T = 0. 

Finite field leads to deviations from the unitary limit 
S a ,b = 7r/2. The magnitude of the deviations is controlled 
by two parameters, B/T a and B/Tb, where T a > Tb are 



the Kondo temperatures corresponding to the two ex- 
change constants. For S = 1 and B = the dot's 
spin is screened completely. If the Zeeman energy is 
small, B <C Tb, then the phase shifts are calculated by a 
st raig htforward application of the Fermi-liquid technique 
P,[l0|, and we find from Eq. (^|) 
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For S > 1 the screening is incomplete: at the fixed point 
da.b = tt/2, and the remaining localized spin is (S — 1). 
The fixed point is reached at B = 0, and the approach 
to it is governed by the ferromagnetic interaction of the 
remaining spin with itinerant electrons fiof . The corre- 
sponding corrections can be calculated with the logarith- 
mic accuracy |ll[] , yielding the asymptote of the conduc- 
tance in the regime B <C Tb, 
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In the weak coupling regime, B 3> T a , one readily obtains 
for arbitrary S ^ (also with logarithmic accuracy) 
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The asymptotes (0), (||) and (^|) clearly demonstrate that 
the dependence of G(B) is non-monotonous with a max- 
imum at B = Tq where Tb < To < T a . This maximum 
can be studied in detail if Xj, <C T a . In this case, the 
intermediate (Tb <C B <C T a ) asymptote for G(B) is 
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and it displays maximum at To = \/T a Tb. The conduc- 
tance is expected to be featureless in the crossover regions 
(B ~ T a or B ~ Tb). The dependence of the conductance 
on the magnetic field is shown schematically in Fig. 0. It 
can be cast into the form 



G = G Fs (B/T a , B/Tb) 



(11) 



The asymptotes of F$ are given in Eqs. (R)-(ff3). 

The dependence of the conductance on temperature 
G{T) is very similar to G(B). The asymptotes for G(T) 
in the leading logarithmic approximation can be obtained 
from Eqs. (p|)-(|l0|) by replacing B — > T, and, in addi- 
tion to it, (S - l) 2 -> S(S - 1), S 2 -> S(S + 1), and 
(5 - 1/2) 2 -» S 2 - 1/4 in Eqs. (|), (§), and &, re- 
spectively. The asymptote for G(T) in the Fermi liquid 
regime (S =1, T <C T b ) follows from Eq. @) after the 
substitution B — > irT. The derivation of these asymp- 
totes is straightforward, but rather lengthy p2^ , since 
the conductance at a finite temperature is not expressed 
via elastic scattering phase shifts [cf . Eq. (0)] . 
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FIG. 1. Magnetic field dependence of the Kondo contribu- 
tion to the conductance through a lateral quantum dot with 
S > 1/2. The temperature dependence of the conductance 
exhibits a similar behavior. 

For the model described by Eq. (|j) conductance G van- 
ishes in the limit J 7 — > 0. Yet, it is well-known that G 
may remain finite even in this limit due to the cotun- 
nelingprocesses, not involving the spin degrees of free- 
dom g. These processes are described by adding to the 
Hamiltonian (Ji|) a potential scattering term Hy- The 
combination of the Kondo and potential scattering yields 
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instead of Eq. (|1 1|) . Here we sketch the derivation of 
Eq. (|l^) for a simplified model in which Hy is diagonal 
in the a — b basis, 



(13) 
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(the general case will be discussed elsewhere Q). Po- 
tential scattering results in finite phase shifts, <5 7S = <5 7 , 
at J 7 = 0; the corresponding clastic cotunneling conduc- 
tance is G e i — Go sin 2 (<5" — 5f), see Eq. In order 
to consider now J 7 ^ 0, one can rewrite the Hamiltonian 
H + Hy in the new basis || , which accounts for the phase 
shifts 8®. In this basis, the Hamiltonian again assumes 
the form of Eq. (||) , but with the modified M exchange 
constants J 7 = J 7 cos 2 6° , which define new Kondo tem- 
peratures T 7 . The phase shifts with respect to the new 
basis are 6 1S — —sS 7 [cf. Eq. (H)]. The phase shifts 5 7 de- 
pend on the magnetic field B via two parameters, B/T a 
and B/Tt>- Returning to the original basis, one finds ||] 
for the phase shifts S 7S = <5 7 — s<5 7 . Substitution of (5 7S 

into Eq. (§) then yields Eq. © with G Q = G - 2G eh 
Microscopic calculation jL2j performed for two specific 
models, the Anderson model and the model of almost 



open dot-lead junctions |13|, gives Go > 0. 

So far, we assumed that both exchange constants J 7 
in the Hamiltonian (Q) differ from zero. Such an assump- 
tion was made in a number of papers, see, e.g., Ref. p4[ . 
However, if a quantum dot is modelled by the Anderson 
impurity then the corresponding effective Hamilto- 
nian (H) inevitably has only one non-zero exchange con- 
stant ]g|8). We will demonstrate now that in fact both 



exchange constants are finite, if Eq. (Q) is derived directly 
from a generic model of a quantum dot. The microscopic 
Hamiltonian of the system consists of the part describing 
the free electrons in the leads: 
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the tunneling part: 
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and the Hamiltonian of the dot: 



(14) 
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The first term in the r.h.s. of Eq. jl^ ) is the single- 
particle part of Hd- The second term, Hi nt , represents 
the electron-electron interaction within the dot. Exten- 
sive numerical renormalization group calculations based 
on such a model were performed in Ref. |l5| , where a non- 
monotonic temperature dependence of the conductance 
was found for a special choice of Hi nt and t an . How- 
ever, the study |L5| did not reveal the two-stage charac- 
ter of the Kondo screening, which shows that the non- 
monotonic behavior is a general feature, rather then a 
result of model assumptions. 

Hamiltonian Hd commutes with the total number of 
electrons in the dot, N = J2 ns 4s^«s, and with its total 
spin, Stot = J2nss' di s (er ss < /2)d ns > (we neglect the weak 
spin-orbit interaction here). The ground state of Hd is 
characterized by the number of electrons N and by spin 
S. In the absence of the magnetic field, the ground state 
is (2S* + l)-fold degenerate and SU(2)-symmetric. This 
can be represented by introduction of the operator of spin 
of the dot in the ground state [cf. Eq. (0)], 



S = VS tot V, 



(17) 



where V is the projector onto the ground state manifold. 

Tunneling Hamiltonian H t mixes the ground state with 
the states having N± 1 electrons on the dot. For brevity 
we assume that the dot is tuned to the middle of the 
Coulomb blockade valley, i.e., its energy is increasing 
by the same amount Eq if the number of electrons is 
changed by ±1. The energy deficit E ~ Ec for the cor- 
responding transitions is high, which allows us to apply 
the second-order perturbation theory in H t (Schrieffer- 
Wolff transformation). The resulting effective Hamilto- 
nian, which is valid at energies less than the gap for intra- 
dot spin excitations A <§; Eq, has the form 
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In the derivation of (|T^) we have neglected the potential 
scattering terms associated with the elastic cotunneling. 

By SU(2) symmetry, the operators VS mn V for any m 
and n should be proportional to S given by Eq. ( |l7| ) : 



"PS V — A S 
This allows us to rewrite Eq. (|l^) as 

H = Hi + ^ Jaa'{s a 'a 1 S), 
act' 

where the matrix of the exchange constants is 
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Note that the factors t an here depend on the proper- 
ties of the tunneling junctions, whereas A mn characterize 
the ground state of the isolated dot, and depend on the 
electron-electron interaction term Hi nt in Eq. (|l6|) . 

We can diagonalize matrix J aa i , and thereby bring the 
Hamiltonian ( |20| ) to the form of Eq. (Q). To do that, we 
perform the unitary transformation ip-yks = J2 a U_y a c a ks, 

The 



Eq. <§. 



where U ia is the same unitary matrix as in 
eigenvalues of the matrix J aa i are J a and Jf, of Eq, 

The matrix A mn of Eq. ([l9|) obviously satisfies the re- 
lations A = At and TrA = 1. For a generic choice of the 
tunneling amplitudes t an these conditions alone fix nei- 
ther the signs nor the magnitudes of the exchange am- 
plitudes J 7 . To make further progress one has to rely 
upon assumptions regarding the form of Hi nt in (16). 
Although model realizations yielding any desired sign of 
J 7 can be easily constructed p2| , here we consider only 
the simplest example, 



H int = E C {N - Af) 2 - E s S 2 t 



(22) 



which nevertheless results in the most interesting case 
J 7 > 0, S > 1. The parameter E$ characterizes the 
strength of the intradot exchange interaction; Es is of 
the order of the single-particle level spacing in the dot. 
In the middle of the Coulomb blockade valley, the dimen- 
sionless gate voltage N is an integer, so the dot in the 
ground state has N = Af electrons. The model (22) is 
fully justified f|,|l6| only for large dots with chaotic single- 
particle states; the relative magnitude of terms which are 
not included in Eq. (p2|), is ~ TV" 1 / 4 <C 1 (we assume here 
ballistic motion of electrons within the dot). 

For the Hamiltonian (|l^), (^2|), occupations of the 
single-particle energy levels are good quantum numbers. 
In the ground state with spin 5, there are 2S singly occu- 
pied energy levels, while all other levels are either empty 
or are occupied by a pair of electrons with opposite spins. 
The matrix A is diagonal, and its only non-zero elements 
are the ones corresponding to the singly occupied levels, 
A m „ = {\/2S)5 mn . Using Eq. (pl|), it is easy to calculate 
the determinant of the matrix J aa i : 
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For S > 1/2 the sum in the r.h.s. here contains terms 
with m ^ n. For a generic choice of tunneling ampli- 
tudes, these terms are positive, and J a Jb > 0. The trace 
of J aa ' , 
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is also positive. Therefore, eigenvalues J a ,b are positive, 
and at S > 1/2 both channels are coupled to the dot via 
antifcrromagnetic exchange interaction. 

In the case of model interaction (^) and the dot spin 
S = 1/2, only one discrete level is singly occulied. Then 
Eqs. (H) and (||) yield J a > and J b = 0, just like 
in the Anderson model. Accounting for corrections || 
to the Hamiltonian ( |22| ) results in finite values of both 
constants, |J&| ~ J a N^^^- 

The two characteristic temperatures, T a and Tf,, should 
be of the same order to allow for an observation of the 
predicted here non-monotonic temperature dependence 
G(T). Because of the exponential dependence of T a ^ 
on the corresponding exchange constants, the two tem- 
peratures are very sensitive to the values of the tunnel- 
ing amplitudes for specific discrete levels, see Eqs. (|2l]), 
(p3|), and (pi|); these values hardly can be predicted in 
advance. However, one can tune the amplitudes by ap- 
plying a weak magnetic field perpendicular to the plane 
of the dot. (It is sufficient to thread about one flux quan- 
tum through the area of the dot p7| ; in the case of GaAs, 
such a field has a negligible Zeeman effect.) The contin- 
uous tuning should allow one to go from a conventional 
manifestation of the Kondo effect at Tt <C T a , to a 
non-monotonic G(T) at T& < T a , to a total suppression 
of the Kondo contribution to the conductance at some 
value of the field, where T b = T a Jl|. 

To conclude, in this paper we considered electron 
transport through lateral quantum dots. We predict that 
if the spin of the dot S > 1/2, then the dependence 
of conductance on temperature and in-plane magnetic 
field is non-monotonic. This work was supported by NSF 
Grant DMR-9731756. 
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